function y = computeTractionLawDerivate(element,etilde) %compute dg/detilde
ef = element.ef;
ft = element.ft;

%exponential law
y = -ft/ef*exp(-etilde/ef);

%linear law
% if etilde < ef
%     y = -ttilde/ef;
% else
%     y = 0;
% end

end
